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Abstract. The extragalactic background of gravitational waves produced by tri- axial 
rotating neutron stars was calculated, under the assumption that the properties of the 
underlying pulsar population are the same of those of the galactic population, recently 
derived by Regimbau & de Freitas Pacheco (2000). For an equatorial ellipticity of e 
= 10~ 6 , the equivalent density parameter due to gravitational waves has a maximum 
<0 amplitude in the range 2x 10~ n -3x 10 -9 , around 0.9-1.5 kHz. The main reasons affecting 
the theoretical predictions are discussed. This background is comparable to that produced 
' by the "ring-down" emission from distorted black holes. The detection possibility of this 
1^ | background by a future generation of gravitational antennas is also examined. 
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1. Introduction 

In the past years, a large number of papers devoted to stochastic backgrounds of gravita- 
tional waves appeared in the literature (see Maggiore 2000 for a recent review). Besides 
processes that took place very shortly after the big-bang, the emission from a large num- 
ber of unresolved sources can produce a stochastic background. Supernovas (Blair et al. 
1997) and distorted black holes (Ferrari et al. 1999a; de Araujo et al. 2000) are examples 
of sources able to generate a shot noise, while a truly continuous background could be 
produced, for instance, by the " r-mode" emission from young and hot neutron stars 
(Owen et al. 1998; Ferrari et al. 1999b). Detection of such backgrounds may probe the 
cosmic star formation rate up to redshifts of z ~ 4-5, the mass range of the progenitors of 
neutron stars and black holes as well as the initial angular momentum of these objects. 
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2 Regimbau.: Gravitational Background 

The contribution of the entire population of rotating neutron stars to the continuous 
galactic background of gravitational waves was considered by different authors (Schutz 
1991; Giazotto et al. 1997; de Freitas Pacheco & Horvath 1997) and, more recently, this 
subject was revisited by Regimbau & de Freitas Pacheco (2000, hereafter RP00). In the 
latter, the "true" population of rotating neutron stars was synthesized by Monte Carlo 
techniques and its contribution to the galactic background of gravitational waves was 
estimated. If the planned sensitivity of the first generation of laser beam interferometers 
is taken into account (in particular that of the French-Italian project VIRGO), then the 
simulations by RP00 suggest that only few objects will contribute to the signal, if the 
mean equatorial ellipticity of neutron stars is of the order of e = 10~ 6 . Upper limits on e 
have been obtained by assuming that the observed spin-down of pulsars is essentially due 
to the emission of gravitational waves. In this case, one obtains e < 10~ 3 for "normal" 
pulsars whereas recycled or rejuvenated pulsars seem to have equatorial deformations 
less than 10~ 8 . Although the galactic population will not produce a truly background, it 
remains to be investigated the integrated contribution of these objects in a fairly large 
volume of the universe. Estimates of this emission are important because it may rival 
with a possible background of cosmological origin. Since there is an upper limit to the 
wave frequency of the pulsar gravitational radiation, a putative cosmic background will 
dominate at low frequencies and, the knowledge of the spectral energy distribution of the 
background produced by discrete sources may help in the choice of the best frequency 
domain to search for a relic emission. 

In the present paper, the integrated gravitational emission of rotating neutron stars to 
the background is calculated under the assumption that the distributions of the rotation 
period and magnetic field derived by RP00 can conveniently be scaled to other galaxies. 
The plan of this paper in the following: in section 2 the model computations are described; 
in section 3 the results are discussed and finally, in section 4 the conclusions are given. 

2. The Model 

The main working hypothesis of our computations concerns the true rotation period 
and magnetic field distributions of pulsars. For galactic objects, these distributions were 
derived by RP00 using Monte Carlo simulations to reproduce the different observed 
distributions of physical parameters, like the period and its first derivative as well as 
distances, when selection effects are taken into account. These simulations permitted to 
establish the parameters of the initial distribution of period and magnetic field. "Faute 
de mieux" , we assumed here that pulsars are born everywhere with rotation periods and 
fields obeying the same distribution laws. 



For a single pulsar, the frequency distribution of the total emitted gravitational energy 
in the source's frame is 
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(1) 



It is worth mentioning that in spite of (dE/dt) in the above equation be the pulsar 
gravitational wave emission rate, the time variation of the frequency is fixed by the 
magnetic dipole emission, responsible for the deceleration of the star. This means that 
angular momentum losses by gravitational waves will never overcome those produced by 
magnetic torques, which is equivalent to say that the average equatorial deformation is 
always less than 10 -3 . Under these conditions, the energy frequency distribution is (the 
fact that the gravitational wave frequency is twice the rotational frequency was already 
taken into account) 

where I is the moment of inertia of the star, r m is the magnetic braking timescale (see 
RPOO), P is the initial period of the pulsar and the other symbols have their usual 
meaning. 

In order to estimate the average ratio {r m /P%), we adopted the following procedure. 
We have performed Monte Carlo simulations in which the distribution probabilities of 
the variables r m and P D are the same as RPOO. The resulting distribution of the quantity 
\og(r m /P^) is given in fig.l, and it can be fitted by a Gaussian with a mean equal to 
<log(r m /P^) > = 12.544. Adopting this value as representative of the whole population, 
the constant K in eq. (2) is K = 9.27x 10 35 (e_6) 2 crg.Hz~ 4 , where we have introduced 
the notation e_ 6 = (e/10~ 6 ). 

The gravitational wave flux at frequency v Q (observer's frame) due to sources localized 

in the redshift shell z,z+dz is 

dv 
dv Q 

where di = (l+z)r is the distance-luminosity r is the proper distance and the observer's 
frequency v a is related to the frequency v at the source by v = (1 + z)v a . The event rate 
inside the shell z,z+dz is 

dR{z) = X p R c (z)^dz (4) 
dz 

In the above equation, R c (z) is the "cosmic" star formation rate, X p = £(m)dm is 
the mass fraction of formed stars in the range 10-40 M Q , supposed to be the mass range 
of the pulsar progenitors, with £(m) being the initial mass function. For a Salpeter's law 
(£(m) oc to -2 - 35 ), \ p = 4.84xl0~ 3 Mq 1 . The element of the comoving volume is 
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dR{z) (3) 
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where H Q is the Hubble parameter and the function E(f^, z) is defined by the equation 

E(Qi,z) = [(1 + zf{\ + zVL m ) - z(2 + z]Q v ] 1/2 (6) 

where Cl m and Q v are respectively the density parameters due to matter (baryonic and 
non-baryonic) and the vacuum, corresponding to a non-zero cosmological constant. The 
equivalent density parameter due to the spatial curvature satisfies Ofe = 1 — Q m — Q v . 
Combining these equations, the expected gravitational wave flux at the frequency v Q 

is 

F Vo KvM—Xe-e) ^ ^.,) dz ( 7 ) 

In the literature is often used an equivalent density parameter due to gravitational waves 
to measure the strength of the background at a given frequency, defined by the equation 

fW - (3^)— (8) 
2.1. Numerical calculations 

In order to evaluate numerically eqs. (7)-(8), it is necessary to specify the cosmic star 
formation rate R c (z) and the parameters of the world model, namely, the values of H OJ 
fl m and fl v . 

Madau & Pozzetti (1999) have reviewed the constraints imposed by the observed 
extragalactic background light on the cosmic star formation rate (CSFR). They con- 
cluded that after an extinction correction of A1500 = 1.2 mag (A2800 = 0.55 mag), a star 
formation rate given by the relation 

- (44 7 + e 3 8 -) M °^- 1 - M ^- 3 W 
fits well all measurements of the UV-continuum and Ha luminosity densities from the 
present epoch up to z = 4. However, according to Hopkins et al. (2001), even when red- 
dening corrections are taken into account, significant discrepancies still remain between 
the CSFR derived from UV-Ha measurements and those derived from far-infrared and 
radio luminosities, which are not affected by dust extinction. Hopkins ct al. (2001) as- 
sumed a reddening correction dependent on the star formation rate and obtained a good 
agreement between the CSFRs derived from different set of measurements. We have fitted 
their results by a function similar to eq.(9), namely, 
1 207p 3836z 

R ^= (39.970 4-^) ^-^0-3 (10) 
Taking into account the uncertainties still present in the derivation of the CSFR, we have 
performed calculations using both rates. 

Recent BOOMERANG and MAXIMA results (dc Bernardis ct al. 2000; Hanany et 
al. 2000) on the power spectra of the cosmic microwave background and observations of 
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distant type la supernovas (Perlmutter et al. 1999; Schmidt et al. 1998), which suggest 
that the expansion of the Universe is accelerating, support a spatially fiat geometry 
and a non-zero cosmological constant. Both set of data are consistent with Q m = 0.30 
(including baryonic and non-baryonic matter) and Cl v — 0.70, which will be adopted in 
our computations. However, no significant differences in our results were observed if a 
model defined by Q m — 1 and Q v = is adopted. The Hubble parameter H Q was taken 
to be equal to 68 km/s/Mpc (Krauss 2001). 

RP00 assume in their simulations that the maximum rotation frequency of a newly 
born pulsar is 2000 Hz, which corresponds to a gravitational wave frequency of 4000 Hz. 
If the upper limit of the integral in eq. (7) is Zmax — 5, then the maximum frequency seen 
by the observer is 660 Hz. For higher frequencies, only near objects will contribute 
to the integral and the upper limit should be replaced by z max — (4000/V o ) — 1, with 
v a in Hz. This parameter affects the resulting spectrum as we shall see below. Thus, 
calculations with a different cutoff were also performed. 

3. Results 

Figure 2 shows the density parameter Qgw as a function of the frequency. Labels Ml 
and HI correspond to star formation rates given respectively by eqs. (9) and (10) and a 
maximum gravitational wave frequency equal to 4000 Hz. The labels M2 and H2 have 
the same meaning but here the maximum gravitational wave frequency cutoff is at 2000 
Hz, corresponding to a minimum rotation period of 1ms. All these curves were calculated 
for an equatorial deformation e=10~ 6 , and we recall that the results scale as e 2 . The nu- 
merical calculations indicate a broad maximum around 1.5 kHz, if the maximum possible 
pulsar rotation frequency is 2 kHz. Decreasing this limit by a half, the spectrum narrows 
and the maximum shifts toward lower frequencies (^0.9 kHz). The amplitude is also 
affected, being reduced by almost one order of magnitude. This happens because accord- 
ing to eq.(8), the amplitude grows as u% but a lower frequency cutoff implies that only 
nearby objects will contribute to the integrated signal and thus reducing the amplitude 
at maximum. 

If the equatorial ellipticity may reach values of the order of 10~ 6 , then the energy 
density of the background generated by pulsars may be comparable and even higher 
than that expected from newly born black holes (Ferrari et al. 1999a; de Araujo et al. 
2000), originated from the collapse of massive stars (M > 40 Mq). For a comparison, the 
spectrum corresponding to the ring-down emission from distorted black holes calculated 
by Ferrari et al. (1999a) is also plotted in fig. (2), appropriately scaled to the Hubble 
parameter here adopted. In the case of distorted black holes, the uncertainties on the 
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estimates of the background energy density rest on the conversion efficiency of the mass 
energy into gravitational waves as well as on the minimum mass of the progenitor. 

Hot and fast rotating newly formed neutron stars may be unstable against the r-modc 
instability. Ferrari et al. (1999b) estimated that if all newly born neutron stars cross the 
"instability window" (see, for instance, Andersson et al. 2000), then the resulting density 
parameter has a maximum amplitude of flow ~ 2 x 10~ 8 in the frequency range 0.5-1.7 
kHz. This signal by far would be the dominant component of the background at these 
frequencies. However, according to the simulations by RP00, only few pulsars are born 
within the instability window, reducing the amplitude of the background due to such a 
mechanism by orders of magnitude. 

Unless the equatorial ellipticity of pulsars be substantially higher than the present 
expectations, the background generated by rotating neutron stars will hardly be detected 
by the present generation of laser beam interferometers and/or resonant detectors, but 
this could be a possibility for future projects presently under consideration, as the Large 
Scale Cryogenic Gravitational Wave Telescope (LCGT), sponsored by the University of 
Tokyo and the European antenna EURO (W. Winkler, private communication). The 
former, with a baseline of 3 km, is expected to have a 100W laser and cooled sapphire 
mirrors among other technological improvements. Therefore, one may expect that ad- 
vanced laser beam interferometers may attain in a near future a sensitivity around 1 



kHz, corresponding to a strain noise \J S n {v) of about 10~ 25 Hz -1 / 2 . On the other hand, 
the best strategy to detect the signal, when the detector output is dominated by the 
noise, which is the present case, is to correlate data from two different gravitational 
antennas and to assume that they have independent noise. One interesting possibility 
would be to correlate the output of such an advanced detector with a resonant mass 
detector located at the same site, having a spherical or truncated icosahedron geometry. 
The advantages of this geometry with respect to a longitudinal bar is that a free elastic 
sphere has five degenerate quadrupole modes, each of which is sensitive to a different 
polarization and wave direction. Moreover, for a given material and resonant frequency, 
a spherical detector has a cross section larger than a cylindrical one. The sensitivity of 
resonant spheres is limited by Brownian motion noise associated with dissipation in the 
antenna and transducer, as well as by the electronic noise from amplifiers. In this case, 
the strain noise at resonance is approximately (Coccia & Fafone 1997) 
AkT 



where k is the Boltzmann constant, T is the sphere temperature, F„ is a dimcnsionless 
coefficient depending on each quadrupole mode (Fi = 2.98, F2 = 1.14, F3 = 0.107), Q n 
is the quality factor of the mode, M s is the mass of the sphere, v s is the velocity of the 
sound and v n is the mode frequency. For practical purposes, let us consider a sphere 
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constituted of the aluminium alloy Al-5056. This material has a sound velocity of 5440 
m.s^ 1 and Coccia et al. (1996) have reported Q values as high as 10 8 for temperatures 
below 100 mK. A sphere with a diameter of 3.5 m (mass of 60.3 tons) has the two 
main frequencies of the quadrupolc modes at 0.8 kHz and 1.5 kHz, covering quite well 
the predicted interval where the maximum amplitude of the pulsar background should 
occur. Assuming a typical temperature of 20 mK, the expected strain noise derived from 



If Av w 20 Hz is the bandwidth of the resonant mass detector and r is the integration 
time, then the expected optimized signal-to-noise S /N of the correlated outputs is (Allen 



For one year integration, one obtains from the equation above S/N w 0.2, indicating that 
new technology detectors may reach in the future the required sensitivity to detect such 
a signal. 

4. Conclusions 

The contribution of rotating neutron stars to the extragalactic background of gravita- 
tional waves was calculated, under the assumption that the parameters characterizing 
the galactic population of pulsars derived by RP00 are the same everywhere. The am- 
plitude of the equivalent density parameter attains a maximum in the frequency interval 
0.9-1.5 kHz and is in the range 10~ n up to 3xl0~ 9 . The amplitude scales as e 2 and, for 
a given equatorial cllipticity, the main uncertainties in the amplitude are essentially due 
to the cosmic star formation rate and to the rotation frequency limit at the pulsar birth, 
which depends on the equation of state of the nuclear matter. For "realistic" equations 
of state, these limits are in the rotation period range 0.5 - 1.0 ms, values adopted in our 
calculations. 

The present estimates indicate that this background, having a duty cycle (measured 
by the product between the typical duration of the signal and the mean birth frequency 
of pulsars) greater than one, may have an energy density comparable to that produced 
by "ring-down" black holes. This emission is unlikely to be detected by the present 
generation of detectors. Correlated advanced detectors may reach a limit of about flaw ~ 
10 — 10 for a flat spectrum (Maggiore 2000), which is not the present case. However, new 
technology detectors, which are presently under consideration, may attain the required 
sensitivity. In particular, taking into account the low cost of a resonant mass detector 
when compared with that of a laser interferometer, the installation in the same site of a 
"sphere" operating near the maximum predict frequency (~ 0.9-1.5 kHz), could be the 
adequate strategy to detect such a signal in the future. 



eq.(ll) is y/Snfa) « 1.5xlQ- 24 Hz" 1 ^. 
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